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Abstract Leading-twist operators have a remarkable 
property that their divergence vanishes in a free the¬ 
ory. Recently it was suggested that this property can be 
used for an alternative technique to calculate anoma¬ 
lous dimensions of leading-twist operators and allows 
one to gain one order in perturbation theory so that, 
i.e., two-loop anomalous dimensions can be calculated 
from one-loop Feynman diagrams, etc. In this work we 
study feasibility of this program on a toy-model exam¬ 
ple of the theory in six dimensions. Our conclusion 
is that this approach is valid, although it does not seem 
to present considerable technical simplifications as com¬ 
pared to the standard technique. It does provide one, 
however, with a very nontrivial check of the calculation 
as the structure of the contributions is very different. 

Keywords conformal invariance • anomalous dimen¬ 
sions 

PACS ll.lO.Hi • 11.25.Db ■ 12.38.Bx 


1 Introduction 

Calculation of anomalous dimensions of composite op¬ 
erators belongs to the standard tasks of any quantum 
field theory calculation. For example, in quantum chro¬ 
modynamics, anomalous dimensions of leading twist 
two operators govern the scaling behavior of quark and 
gluon distributions in hadrons and have to be known 
with high precision. Nowadays the anomalous dimen¬ 
sions are known at three-loops, see Refs. [ini and refer¬ 
ences therein. Beyond the two-loop approximation such 
calculations are feasible only with the help of the ad¬ 
vanced methods of computer algebra. Since the calcula- 
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tions are fully automated finding errors becomes highly 
nontrivial task and any approach which can provide a 
check of the final results is very helpful. 

Given that the theory depends not only on the cou¬ 
pling constant but also some other parameters such 
as the dimension of an internal symmetry group, one 
can organize expansion over these parameters. The best 
known example of this kind is the I /TV expansion, see 
Ref. for a review. Agreement between the results ob¬ 
tained in the perturbative and \/N expansions serves as 
a powerful test for the validity of calculations. However, 
the calculations in the 1/N expansion are much harder 
than perturbative calculations: only two RG functions - 
indices of the basic fields in the nonlinear tr—model and 
Gross - Neveu model - are available at 1/N^ order 
01 ■ In QGD the calculations rarely go beyond the lead¬ 
ing order in 1/Nf (where Nf is the number of flavors). 
At leading order in 1 /Nf the anomalous dimensions of 
twist two operators were calculated in Refs. 0 but 
extension of these results to the next order is hardly 
possible. 

A new approach for calculating the anomalous di¬ 
mensions of leading twist operators was proposed in 
Refs. l3, 111- It is still a perturbative approach, how¬ 
ever, the contributing diagrams are completely different 
from those in the standard technique. The approach is 
based on a remarkable property of leading twist opera¬ 
tors: namely, a divergence of such an operator, , 

vanishes in a free theory 121 


( 1 ) 

In the interacting theory the r.h.s. of Eq. o is non¬ 
zero but is proportional to the coupling constant. This 
identity allows to extract the loop contribution to 
the anomalous dimension of the operator from 

i — 1 loop diagrams only. In particular the one loop 
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anomalous dimensions of leading twist operators do not 
require calculation of loop integrals at all 


11, 2C 


The method developed in Ref. [ll| is adjusted to lo¬ 
cal operators and relies heavily u pon the so-called ” con¬ 
formal scheme” renormalization l3“ 15|. In our opinion 
it is more convenient to stay within the standard MS 
scheme and the formalism non-local (light-ray) opera¬ 
tors technique. This technique proves to be more effec¬ 
tive and flexible as we will demonstrate on the example 
of calculation of two - loop anomalous dimensions in 
the su{n) symmetric (p^ model jl^ . 

The paper is organized as follows. In section [5] we 
introduce the model and fix notations. In section |3] we 
recall the light-ray operator technique. Section H] is de¬ 
voted to calculation of the divergence of conformal op¬ 
erator. The details of calculation of two - loop corre¬ 
lators are presented in section Our conclusions are 
in section [T] In Appendices we explain some technical 
issues and details of the derivation. 


2 Generalities 


The su{n) symmetric (p^ model is a scalar field theory 
in d = 6 — 2 e = 2 /x dimension with an action 


S{p) = j df'x 
where a = 1 ,..., — 1 , 


+ , ( 2 ) 
z o 


= 2trt“{t'’t“}, 


(3) 


are the generators of the su{n) algebra normalized 
in the conventional way, trt“t^ = 1/2. The theory is 
multiplicatively renormalizable 


SR{p,g) = S{pQ,go) , 


(4) 


where 


ipQ = ip and go = M^Zgg . 

Two-loop expressions for the renormalization constants 
Zi = and Z^ = ZgZ^ can be found in Ref. [l^. The 
/3—function of the charge u = g^/( 47 r)^ and the field 
anomalous dimension 7 ,^ are 

j5{u) = —2eu — + 0{u^), 

2n 

, , - 4 / - 100\ 3 

7^M = »^^(i+»^^)+0(. ). (5) 

At the critical point u*, (3(u^:) = 0, 

u* = 4ne/(20 — n^)-I-O(e^) ( 6 ) 


theory enjoys the scale and conformal invariance. 0 
It is well known that in a conformal theory the form 
of two-point correlation function of conformal opera¬ 
tors is fixed up to normalization. In particular, the cor¬ 
relator of the conformal traceless symmetric operators 
has the form 


(pf\x)0f{y)) = 6,,, 6 a, A, 


{{x - yY)'^i 


(7) 


Here j {j') is the spin of the operator. The vectors n 
and n are two light-like vectors, in? = fi^ = 0 , and 
0^^\x) {0'^^\y)) is a contraction of an operator with 
the vector n (n), for instance 


0l"^\x) (8) 

Aj and Aj/ are the scaling dimensions of the operators, 
Cj is a normalization constant and 


In,n{x) = (nn) - 


2{nx){nx) 


(9) 


The scaling dimension of the operator is given by the 
sum of canonical and anomalous dimensions at a critical 
point, Aj = A^^'^ + 7 j, where ■jj = ^j(u*)- For the 
leading twist operators of spin j, = 2fi — 2 + j. 

The anomalous dimension of an operator in MS 
scheme, 'Jjiu), is a function of a coupling constant only. 
It can be restored from its critical value, 7 ^ = 7 j(u*), 
provided that the latter is known as a function of e. 

The correlation function for the divergence of the 
conformal operator 

doY^{x) ( 10 ) 

can be obtained from Eq. ([7]). If x is chosen in the trans¬ 
verse plane, (xn) = [xn) = 0 , the ratio of the two cor¬ 
relation functions 


'Tj{u?) 


x^{nn) 


(dO^J"\x)dOf\Q)^ 

iof\x)of\Q)) 


( 11 ) 


is a function of the anomalous dimension 7 ^- only 


T)- = 2 / 7 j 


{2p. - i + j){p - \ + j) 


p — 2+ j 


2jj 


( 12 ) 


This ratio, in full agreement with the result of Ref. [l^ . 
vanishes provided that 7 ^ = 0 . 

^Formally, a nontrivial critical point for d < 6 only exists for 
n = 3, 4. However, staying within perturbation theory one 
can consider n as a continuous parameter. In this sense all 
further results hold for arbitrary n. 
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The perturbation series for 7 ^ and Tj 
( 1 ) , 2 ( 2 ) , 

Ij = M* 7} + + ■ • ■ , 

7-=^, + , (13) 

are related to each other. For later convenience we chose 
the normalization factor Kj as follows 


2j-0 '+m + 3) 
j +1 


(14) 


Substitutiii^ tliG SGri 0 S into m one finds the 
following relations between the expansion coefficients: 


7 . 


( 1 ) 


( 1 ) 


„( 2 ) 


= T, 


_ 7 ^)^) _ 


+ 2)(j+3)^ ^ ^ 


+ 


(j 

2f + 5j + l 


(j + l)(j+2)(j+3) 


71) 71(1) 


(15) 


and so on. Here = (n^ — 20)/4 ?t,. 

Since the divergence of the conformal operator is 
proportional to the coupling constant, dOj^^ ~ 0{g), 
the ratio Tj contains a ” kinematical” factor u . 
Thus, in order to determine Tj and, hence, the anoma¬ 
lous dimension 77 , with 0{u^) accuracy the correspond¬ 
ing correlation functions have to be calculated at one 
order in u less. In Ref. [lll| one loop anomalous dimen¬ 
sions were reproduced by this method in and = 4 
SUSY models. Going beyond the leading order requires 
an effective technique for calculation of the two-point 
correlation functions otherwise one gains nothing in 
comparison with the standard approach. 


Here [Okm\ is the renormalized (in MS scheme) local 
monomial Okm = d^ip°‘{x)d'^g}°'{x)lk\m\ and 9+ = 
{ndx)- The sum in Eq. (flTll can be replaced by action 
of some integral operator on the bare operator 


[0(x; z)] = ZO{x; z). 


(18) 


Here we introduced a shorthand notation 2 ; = {zi, Z 2 
The integral operator Z can be written in the form [1 


Zf{z) = J dadpZ{a,f3)f{z‘^^,z^,). 


(19) 


The renormalization kernel Z{a,f3) is given by a series 
in 1 /e and the coupling u. The light-ray operator (fTTll 
satisfies the RG equation 

(^MOm + Piu)du + M{u)j [0{x; 2 )] = 0, (20) 

where the evolution kernel H is an integral operator 

m{u) = - + 27 ^ ( 21 ) 

which encodes all information on the anomalous dimen¬ 
sion matrices for local operators. 

At the critical point u = local operators can be 
classified according to the representations of the con¬ 
formal group. An operator with the lowest scaling di¬ 
mension in the representation is called a conformal op¬ 
erator. The leading twist operator H, Oj, is uniquely 
determined by its scaling dimension Aj. The expansion 
of the light-ray operator over conformal operators 
and their descendants reads [l^ 


3 Light-ray vs local operators 


The first task is to find a convenient description for 
local operators. As we will argue the light-ray opera¬ 
tor technique 17| is a most suitable one. The light-ray 
operator 


[0{x; zi,Z 2 )] = + zin)(p°-{x + Z 2 n)] 


(16) 


is defined as the generating function for the renormal¬ 
ized local operators 

[ 0 {x-,Zi,Z 2 )\ =^z\z'^[Okm]{x) 

k,m 

= Y. zlz^Zf^r'Ok'm'ix). (17) 

k,m,k' ,m' 

^In order to make presentation more transparent we will con¬ 
sider the su{n) scalar operator. The operators of other sym¬ 
metry properties can be easily included into consideration, 
see Sect|^ 


[0{x- z)] = Y '^ 3 k{zi,Z 2 ) dXO,{x ), ( 22 ) 

jk 


where the coefficients ’F,fe(zi, Z 2 ) are homogeneous poly¬ 
nomials of degree j -I- fc in zi, Z 2 • These polynomials are 
eigenfunctions of the evolution kernel ]HI(m») 

lHI(M,)!F7fe(z) = 77 'I'jkiz ), (23) 

with the corresponding eigenvalues being the anoma¬ 
lous dimensions, 77 = Since the theory enjoys 

conformal invariance at the critical point u = the 
evolution kernel commutes with three generators of the 
collinear subgroup of conformal group, 0 


[^±.o,H(u,)] =0. 


(24) 


®Note, that the operator Oj vanishes identically for odd j. 

^ Other generators acts on the operators in question trivially. 
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The generators, however, deviate from their canonical 
form (see Appendix [Appendix A[ ) 

= -a.i - d,, , 

Sq ^ = zidz^ + Z2dz2 + 2 , 

= zfdzi + + 2(zi + Z 2 ) (25) 


The integration goes over the unit disks |zfe| < 1, k = 
1,2. The generator is a self-adjoint operator with 
respect to this scalar product and (5’^°^)'!' = —We 
want to find a deformation of the scalar product m 
that keeps these relations for the complete generators. 
So = si and Sl_ = —S-. Let us look for the solution in 
the form 


due to quantum corrections, 

5±,o = 4°0 + ^^±,0- (26) 


Two of the generators are known to all orders, 

AS- = 0, ASo = -€+ iH(n,), (27) 


while corrections to the generator of special conformal 
transformations can be calculated order by order in per¬ 
turbation theory 1^ . The leading correction is 


AS+ — {zi + Z 2 ) ( —e+—l 


)]+0{e^). 


(28) 


It follows from Eqs. (EH), EH) that the operators S± 
act as raising (lowering) operators on the set of eigen¬ 
functions 'Pjk, 


S±<I'jk ^ tfbfezti • (29) 

In turn Sq counts conformal spin of the operator Ojk-, 

So^jk — jjk^jk — 2 i^jk T *^jfc) ^jk : (^6) 

where Ajk = Aj + k and Sjk are the scaling dimen¬ 
sion and spin of the operator, respectively. One derives 
immediately from EH that the polynomial ^jk=o ac¬ 
companying the conformal operator Oj in the expan¬ 
sion (l2^ is a simple power, 'l'jk=o{zi, Z 2 ) ^ (zi — ^2)-^ 
and all other eigenfunctions have the form ^jk{zi, Z 2 ) 
5^(zi-Z2)L 


3.1 Scalar product 


Eq. (l23ll can be considered as a standard quantum me¬ 
chanical problem for the Hamiltonian H. In order to 
make this analogy complete one needs to introduce a 
scalar product on the space of the eigenfunctions. Clearly, 
such a scalar product has to be adjusted to the sym¬ 
metries of the problem. At leading order the answer 
is given by the standard sl(2) invariant scalar prod- 

30 


uct 


(■*/','(')o = ;^ jj<fzi(fz 2 {'ll}{zi,Z 2 ))^(l){zi,Z 2 ). (31) 


bfc|<i 


(^/’,<?i’)ro = (V',ro<(')o , w = 1 ++ ... , (32) 

where is a self-adjoint operator with respect to the 
scalar product ED). The conjugation conditions for the 
generators imply 

ZIS'W = , (33) 


The one loop corrections to the generators involve the 
kernel (EI(u) = which is given by the 

following expression 

^(1) ^ 2(7^ - A«77+). (34) 

Here As is a color factor. As = {'n? — 4)1 n and 


'H+iP{z)= f da [ d^V'( 2 ? 2 >^ 2 i) • 
^0 ^0 


(35) 


The operator is completely determined by Eqs. ED- 
The explicit expression for and details of the deriva¬ 
tion can be found in [Appendix B| [^. 

Since the eigenfunctions Wjk are mutually orthogo¬ 
nal w.r.t. the scalar product EDi on® can represent the 
conformal operator as the scalar product of the coeffi¬ 
cient function with the light-ray operator 


= {z{^,[0{x,z)])zzz 


(36) 


This representation for the conformal operator is the 
most convenient one for further analysis. We demon¬ 
strate it on the following example. The conformal op¬ 
erator is usually defined as an operator which vanishes 
under special conformal transformations, Kn = K ■ fi, 
(5if_(!Ij(0) = 0. This property becomes transparent in 
the representation (IMl) if one takes into account that 


5KAO{z)]=2{nn)S+[0{z)], 

(we put here [0(z)] = [0{x = 0;z)]) and use that the 
generators S+ and S- are conjugate to each other w.r.t. 
the scalar product ED, S'! = -5_. 

®It turns out that the corrections due to cancel at O(e^) 
order in the ratio of the correlators CD. So that we do not 
need this explicit expression for the present purposes. 










5 


4 Divergence of conformal operator 


In order to construct the divergence of the conformal 
operator d0^p\ Eq. (ITOll . we calculate first the diver¬ 
gence of the light-ray operator da 




(37) 


Taking the n—derivative one cannot, however, keep v? = 
0 any longer and has to take into account terms linear in 
v? which account for the trace subtraction in [0{x', z)]. 
Taking the corres pon ding modification into account, see 
e.g. Refs. 11, 17, 211, one gets for (11171) 

[dO{x\ z)] = ^Zip°-{x + zin)ip°'{x + Z 2 n) (38) 
where is a differential operator 




(39) 


which commutes with the renormalization factor Z and 
acts on the fields directly. After a simple algebra one 
gets 


[dOix;z)] = - l)-^zi^S^^^dlO{x-,z) 

— L^i d^ip°‘{x + zin)(p°‘{x + Z 2 n) 

- ip‘"{x + zin)d‘^ip°'{x + Z2n) > , (40) 


where — e, — e{zi + Z 2 ) and 

= '1^22-^21 ~ ^■^21 , Ti2^ = ^212^12 ~ f - Zl2 ■ (41) 

Using equations of motion (EOM) one can replace in 
this expression 

9V“(^) ^ id“''>'’(x)v5“(x). (42) 

We want to stress here that Eq. dini) holds for arbitrary 
coupling u but not only at the critical value. Since the 
l.h.s. of Eq. (HU)) is a finite (renormalized) operator the 
r.h.s. can be expressed in terms of renormalized oper¬ 
ators with finite coefficients. These operators can be 
chosen as: the two-particle operator Oi = d^O(x; z) 
and three particle operator 


where w = {wi, UI 2 , W 3 }. The operators 0\ and O 2 mix 
under renormalization. The mixing matrix (integral op¬ 
erator acting on fields variables) has an lower triangular 
form 

[Ok] = ZkmOm . (44) 

Here Zn = Z is the renormalization constant of the 
light-ray operator, Eq. (dll), Z 12 = 0 , Z 21 = 0{u^) and 
the element Z 22 is given, at the one loop order, by the 
sum of two - particles kernels 

^11 = 1 -^ A, 77+ + 0 (^ 2 ), 

^22 = 1 + ^ ’ (45) 

where , \d are color factors 

As = (n^-4)/n, Ad = (n^ - 12)/n. (46) 

The kernel Hff. is defined by Eq. (IHK|) and Jiff, has the 
form 

77 ^ 2 /( 21 , ^ 2 ) = / da aa f{z ^2 0 12 ) ■ (47) 

The subscripts ik show the arguments the kernel acts 
on. 

Using these results we can rewrite dMl) as follows 

[dOix-,z)] = A,[0,(x;z)]. (48) 

The operators Ak have the following form: 

Ai = - M-M 2 Z 2 i^n , 

( 72 ^ 2^^2 + Ui 5i) Z 22 I, (49) 

where the operators Si, S 2 map functions of three vari¬ 
ables to functions of two variables 

[Slf]{zi,Z 2 ) = f{zi,Zi,Z 2 ) , 

[S2f]{zi,Z2) = f{zi,Z2,Z2) ■ (50) 

At one loop the operators Ak take the form 

Ai = S+(u) + m(^As77+ - 7 ^ ^)('Zi + ^ 2 ) + O(u^), 

A 2 = ——M^^Li2S2 + T2i5'i -I- — uAs77+^zi 2 >S'i 2 

- ^^12^12 (a« 77^3 - Ad77+3) + 0(^2)) . (51) 

Here <S'+(ti) is given by the expressions ( 1211 ), (EZl) for 
arbitrary u, ^ u, 

Q _ Q Q T _ T (^'^0) 

O 12 — — O 2 , iJkm — J^km 


O 2 = 0^^\x-,w) 

= -|- win)(p’^{x + W2n)ip^{x + wsn ), (43) 


and deriving dSU we made use of the symmetry of the 
three particle operator Q^'^\x]Wi,W 2 ,w^) under per¬ 
mutation of w variables. 
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Let us stress again that the operators Ak do not 
contain singular terms for arbitrary u. Using one loop 
expressions for Z factors, Eq. gS]), it can be checked 
that all pole terms cancel at order 0{u). In particular 

■^ 3-^1 ^-^11 ^ Li 2 * S '2 + L 21 * S ' i ^ ^ 22 ^ = 

= L 12 S 2 + L2i<5'i + 0{u^). (52) 

Starting from the representation (Ei) for the con¬ 
formal operator, we get for its divergence 

dOj{x) = (zf 2 , [dO{x-,z)])^ . (53) 

Making use of Eqs. (H51) - (1511) one finds that the diver¬ 
gence dOj(x) is given by the sum of two-particle and 
three-particle (renormalized) operators 

= + (54) 

The prefactor on the r.h.s of Eq. (15^ is the eigenvalue 

of the operator (S'q'^ — 1)“^ on the function z{ 2 - The 
( 2 ) 

two-particle term TVj has the form 




Fig. 1 The leading order diagrams for the correlator of two 
conformal operators, (0)). 


of Eq. (15^ , results in the following relation for the cor¬ 
relators 

{nf {x)nf (0)) = - {nf {x)nf (o)). (eo) 

It can be shown that in the correlator {dOj{x)dOj{Qi)) 
one can replace dllD by a simpler expression 

dO,{x) = (nf [x) + 7^f (x)) . (61) 

The omitted terms 


Tlf\x) = {zi 2 , {zi + Z 2 )[Oi(x;z)]} + O(m^). (55) 

We recall that jj = 7 j(u*) and Oi{x;z) = d^O{x-,z). 
Let us note that the term ^ S+ in the expression for 
Ai, Eq. (EH), vanishes inside the scalar product since 
(z^ 2 ) ■ ■ ■) = ~i‘S'-zl 2 ,...) = 0. In turn, the expres¬ 

sion for the three-particle contribution can be written 
as follows 

nf{x) = nf’°\x) + nf’^\x) + oie ^), (56) 

where 

nf’°\x) = -M\z{ 2 , L 2 iS^p 2 {x-z)])^ , 

Tlf^\x) = -M^{z{2,z,2S^X,[02{x-,z)\)^ (57) 

and the operator Xj has the form 

Xj =e-^^+ 7j/ 2 - u* {XsTiis - Xd'Hts) ■ (58) 

This expression follows immediately from m if one 
takes into account that spin j is even and A 2 is sym¬ 
metric under interchange zi Z 2 - 

It is clear from (I55|l that the expansion of two parti- 
( 2 ) 

cle term 77.) ^ over conformal operators does not contain 
the operator of spin j, 

^ d‘^ ■ (59) 

Vm^O / 


“ JTT 

give rise to the correction of order O(e^). In order to ver¬ 
ify this it is sufficient to notice that AX can be rewrit¬ 
ten in the form {z{ 2 t F S-Zi2[02{x] z)]) where F is 
some operator whose explicit expression is not relevant. 
Inside the correlator the generator S- acts, finally, on 
the function 21)2 nullifying it. 

Thus in order to find the anomalous dimension 7^ 
at order O(e^) one has to calculate the three correlators 

at one loop order. We will do it in the next section. 

Finally, we note that for j = 2 the r.h.s. of Eq. dSH) 
has to vanish identically since the operator is, up 
to EOM terms, the energy-momentum tensor.The two 
particle term 77) is proportional to the anomalous di¬ 
mension 7j and therefore vanishes for j = 2. In order to 
(3) 

check it for 77) , it is sufficient to take into account that 
only the linear term in the expansion of three particle 
operator 

02{x;w) ~ (ici -I-UI2 +W3) ■ d°-^‘'d+ip°'{x)(p^{x)ip‘'{x) 

= • l) d“''=a+(/?“(x)(/J^(x)(/?=(x) 


Taking into account that {Oj{x)Ok{0)) = 0 for fc < j contributes to (I5S1) for j = 2. After simple algebra one 

one derives that {TZj^\x) dOj{0)) = 0. This, in virtue finds that 77^=2 = 0{e‘^). 
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Fig. 2 The LO diagrams for the correlator of divergence of 
conformal operators, (90^™^ 


5 Correlators 

5.1 LO correlators 


In order to give a glimpse of the technique we start with 
calculation of the necessary correlators at leading order. 
The correlator is given by the sum 

of two diagrams, shown schematically in Fig. [TJ They 
are given by the product of the propagators and give 
rise to identical contributions to the correlation func¬ 
tion. Assuming that x is chosen in a transverse plane, 
(x,n) = {x,n) = 0, one represents the propagator as 

D{x + zin — Win) = Zl(a;)(l — ziwir)~^^~^'^ , (62) 

where r = 2{nn)/x'^ and 

D{x) = r{n-l)/{47T>^{x^)^-^). (63) 


At leading order one replaces /j, — 1 i-d- 2 so that the 
second factor in (l62ll is nothing else as the reproducing 
kernel, ICs=i{zi,wir), corresponding to the spin s = 1, 
see Eq. (IA.5I) . Therefore starting from Eq. (1551) one gets 
for the correlator 

2 

(o]”^(a;)C>j"^(0)) = 2^D'^{x){zl2\ JCi{zk,Wkr)\wl2) 

k^l 

= 2^D^{x)r^ Wwl^Wli, (64) 

where ^ — 1 is the isotopic factor, the scalar prod¬ 

ucts correspond to the conformal spin s = 1 and we 
take into account the property of the reproducing ker¬ 
nel (IA.6I) . The norm of 111(2 is given by the following 
expression 



2 

S1S2 


.| TT r(2sk) T(2j-b 2(si -b S2) — 1) 
" 1 ^ r{j + 2sk) r{j + 2(si -b S 2 ) — 1) 


The diagrams for the correlator (90j”^(a;)90j"^(0)) 
are shown in Fig. [2] On the leftmost diagram the points 
zin and wfi are connected by two propagators 

D‘^{x + zin — Wifi) = D^{x){l — ziwir)~‘^^^~^^ 

= D^{x))Cs=2izi,wir). (65) 


Since 90^"^ ~ (2(21 .^2i[0(x, zi, zi, Z2)])ii, see Eq. dST]), 
the z— scalar product has the form 

{z{2,L2ilC2izi,wir)ICi{z2,W2r))ii . (66) 

The spins of the reproducing kernels and spins of the 
scalar product are in discord with each other. However, 
the operator L21 = d 2 Z 2 i removes this mismatch. It 
intertwines the representations, 

L 21 D+ (g) D+ = £)+ (g) D+L 21 

and it can be easily shown, see e.g. Ref. [^, that 

(z(2, L2l<P{zi, Z2))ll = -aj(z(^\ ^>(Z1, Z2))21 , (67) 

where 

aj = (J+ l)lk(2ll?i/IN^2-'ll^i = + + (68) 

Thus the scalar product (15511 takes the form 

-aj{z{2^,IC2{zi,wir)ICi{z2,W2r))2i = -aj{rwi2)^~^. 

(69) 

Restoring all color and symmetry factors one gets for 
the first diagram 

- ^^>^sg^D^ix)aj/{j + l)^r^-^ {L2 iw{ 2^ ,w{2)ii = 

= ^i>^sg^D^{x)aj/{j + l)r-^"^ 11^(21 111 • (70) 

The calculation of the second diagram goes along 
the same line. One can combine propagators attached 
to the point zi using the Feynman’s trick to get 


(z(2, L2ilC2izi,wir)IC2(zi,W2r)ICi{z2,wir))ii = 

_ 6(-l)%j-(rwi2)^~^ , . 

(j + l)(j+2) ■ ^ ^ 

Finally, taking into account that second diagram en¬ 
ters with symmetry factor 2 one gets, in full agreement 
with (1TT1) . 

7) (u*) = -b 0{ul ), (72) 

where the one-loop anomalous dimension is, see Eq. (IMl) . 


(1)^ 1 (j - 2)(j + 5) 

O'. «6(j + 1)(j + 2)- 


(73) 


Thus the diagrams are easily calculated provided that 
the spins of the reproducing kernels match that of the 
scalar product. We will show that this scheme can be 
extended to loop diagrams as well. 
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Fig. 3 


NLO correction to the correlator of conformal operators 


{Oj^\x)Oj^\o)). The parameter 5 


e/2. 


Wl 


w 


P 

12 


5.2 NLO correlators 


First, it is easy to see that the corrections due to mod¬ 
ification of the scalar product cancel out in the ratio of 
correlators. Indeed, these corrections only influence the 
norm 


^121111 


I—>■ kal 


= (wi 




entering the tree level expressions, Eqs. (EH), dZQl) , etc., 
which cancels out in the ratio of correlators, 7j, irre¬ 
spectively of the explicit form of w. This cancellation 
is expected. Indeed, the problem can be reformulated 
as a standard quantum mechanical problem for a cer¬ 
tain Hamiltonian. A modification of the scalar product 
produces corrections to the eigenstates (conformal op¬ 
erator). However, the energy shift at the leading order, 
is not sensitive to such correc¬ 
tions. 

The calculation of is a bit more in¬ 

volved but straightforward. The corresponding contri¬ 
bution to reads 




(2^3) 


1/ (i)',2 (j ~ -I-9) 

2^ ^' ^ (j + 1)(j + 2)(j + 3)(j + 5)- 


(74) 


It is convenient to split one loop corrections to the cor¬ 
relators {OjOj) and into two groups 

— Self-energy insertions to the propagators. 

— All other loop diagrams 

Taking into account the self-energy correction to the 
propagators is equivalent to the calculation of the tree- 
level diagrams with the exact (critical) propagator 

Dc{x) = , (75) 


where s = = 1 —(e—7,^)/2, the subscripts indicate 

the conformal spins of the z and w scalar products, re¬ 
spectively. The reproducing kernels in dZZD correspond 
to spin s which does not match the spins of the scalar 
product. However, it is easy to see that, due to symme¬ 
try, the scalar product with modified spins 

(2:(2|/Cs(zi,ri;i)/Cs(z2,W2)m2)(4.4).(si.si) , (78) 

where = 1 ± <5 is equal to that in (1771) up to terms of 
order . Therefore, choosing si = s, (S = (e — ^^)/2), 
and evaluating the ic—product in (17^ one gets for (1771) 

Dl{x)r^\\zUls^^^.^+0{ul). (79) 

The calculation of the diagrams in Fig. [7] goes along 
the same lines. Finally, the contribution to the ratio 
of correlators 7j from the leading order diagrams and 
self-energy diagrams can be written in the form (up to 
O(e^) terms) 


q-SE 
' 3 





^ _ 2T(4sl)(T(j-l + 2sl) \ 

"" V r(2sl)T(j-l + 4sl) y ■ 


(80) 


Expanding (18011 we find for the corresponding contribu- 
( 2 ) 

tion to the coefficient 




2<5 


’2i-|-2 


“ Sj+3 — Si + - 


2 4/2 -I- 14j -I- 9 


3 (j + 2)(j + 3) 


where A^ = ^ — 1 + is the critical dimension of the 
basic field and the residue A(u*) is 

AM = ’ (^6) 

V 36 y 

where M = Eor the first diagram in Fig. [T] 

one gets 

Dl{x) {zi2\Kis{zi,wi)lCs{z2,W2)\w{2){ii),(ii), (77) 


4A,<5 

(J + 1)(J + 2) 


■ (2j + l)(4j + 7) ' 

. 3(/ + l){j + 2) _ 


, (81) 


where Sj = X)fc=i 1/^ 

^ = - -I- ■ (82) 


We recall that is the one loop anomalous dimen¬ 
sion, Eq. dZH. 
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5.3 Loop diagrams 

All loop diagrams can be calculated quite easily with 
the help of several simple tricks. Several examples are 
given below. Let us start with a correction to the cor¬ 
relator of two conformal operators. The corresponding 
contribution has the form 



C(e) = 2 {z{ 2 \H{x; z, w)m 2 )(ii),(ii) > (83) 


where the kernel H(x; z; w) is given by the left diagram 
shown in Fig. [3] with the parameter d 0. We need to 
find C{e) up to terms 0(e°) , (7(e) = i(co-l-eci-|-...). To 
this end we proceed as follows. We modify the indices 
as shown in Fig. [3] This modification does not change 
the pole structure (see discussion in Ref. [26|) and, due 
to the symmetry (7(e, <5) = C{e, —S), one concludes that 


^) — ~(co + -|- C2(5^ -I- ...). 


(84) 


The choice d = e/2 results in the uniqueness of the 
upper integration vertex and at the same time does not 
affect first two terms in (IMl) . Using the star-triangle 
relation for the upper vertex one gets the diagram B 
in Fig. [31 Using Feynman formula for the left (right) 
propagators attached to the integration vertex one can 
perform the last integral that results in the diagram C 
in Fig. [3] 

This diagram, up to a a;—dependent factor, has the 
form 


ICi{zi,wi) / dadj5{aj5) \aPy /C|_35(0j*2,K 2 ) • 


(85) 


On the next step we want to get rid of the paramet¬ 
ric integrals. To this end we use the properties of the 
reproducing kernel (IA.6I) and represent 


K,s-5[Zi2.w(y) = 

J d^e J d"eM.(c)-*Cs-6(c',o^.(e,<) 

= J d^^fis+s{0’^s{z^2,0ICs{^,W^12)+O{6^), ( 86 ) 


Fig. 4 The “sl(2)” diagram: an arrow line from w to z with 
index a stands for the propagator (1 — zw)~°‘. The indices 
have the following values: a = 2 — 3e/2, /3 = 1 — e/2 and 
7=1. The black circle denote an integration vertex with the 
sl{2) invariant measure fis+s, s -|- <5 = 3/2 — e/4. 


scalar products, S{a, b, si, S 2 |a', b', s'^, S 2 ). We need this 
function up to 0{e) terms ior a = a' = a, b = b' = fd 
and Si = s' = 1. As was explained in the previous 
section the integral with measure fj-s can be evaluated 
easily provided that the sum of the indices of propa¬ 
gators coming from this vertex is equal to 2s. Taking 
this into account one finds that the scalar product with 
shifted indices 


S{a -h d, /3 -h (5,1 -h 2(5,1 -h 4d|a — 6, /d — 6,1 — 26,1 — 46) 


can be straightforwardly calculated to 


S{e) 


U(j+3-e) IUi2lli+iM+. 


(87) 


and in the same time it differs from the scalar product 
in question by terms of order 0{e^) only. 

Restoring all factors one gets for (7(e) 

(7(e) = u* As(n^ — 1)D‘^{x){x^M^Y ^ ^ S'(e). (88) 

e 

Since the ratio (HU does not depend on x, it is conve- 
nient to put x M =1. Finally, subtracting the coun¬ 
terterms 


AC(e) = - V A.(u^ - l)D^x)r^ 


one obtains 

C(e) + AC(e) ^ 2(52 ,+2 - 5,+ 1 ) 

{Oyx)OM), ^ (j + l)(j+2) 


(89) 


where s = 3/2—2(5 and we used the relation (IA.7|) . Using 
this expression one can carry out the integrals over a, /3 
in (|85l) . The resulting expression for H(x] z, w) (up to a 
prefactor) takes the form of “si(2)” diagram shown in 
Fig. m 

Now we need to calculate the scalar product (l83|) 
with the kernel H(x;z,w) given by the diagram on 
Fig. m The scalar product is a function of indices of 
left (right) propagators and conformal spins of 2 ; (w) 


where ellipses stand for higher order terms. The corre¬ 
sponding contribution from (1891) to the coefficient 
in the ratio of the correlators, see Eq. (HU), reads 


rpio) ^ 2(52j-s2 - 5,-+i) 

^ (j + l)(j + 2) 


(90) 


The NLO diagrams contributing to the correlator 
{dOj{x)dOj{0)) are shown in Fig. |S] (type A) and Fig. [5] 
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W2 


W\ 


Fig. 5 One loop correction to the correlator of the divergence of conformal operators {x)dO^^'' {Q)). 


(type B). These diagrams have different color factors: 
As Ad for the A - diagrams, and Af for the B - diagrams. 

All diagrams of the A - type have the diagram we 
have just now discussed as a subgraph. The only differ¬ 
ence is that in order to kill terms linear in <5 one has to 
consider an average of the diagrams, {D{5)+D{—5))/2. 
For 5 = e/2 each of the diagrams, D{±d), can be sim¬ 
plified with the help of the star-triangle relation and 
rewritten in the form of “sl{2)” diagrams. These dia¬ 
grams in turn can be calculated up to O(e^) terms in 
a manner described above. So we skip all details and 
present the result for each diagram in [Appendix C[ 

All diagrams of B-type shown schematically in Fig. [5] 
contain two 2 —1 subgraphs. The diagrams depicted 
on the right panel are finite while those on the left panel 
are divergent. The calculation of these diagrams does 
not present any problem so that we give only final re¬ 
sults in [Appendix C| 

Finally, it follows from Eq. dSll) that the sum of 
counterterm diagrams to the diagrams in Figs. [5] and [B] 
can be written in the form 

2{Z3Z^^Zj - 1) {dO^"\x)dOl"\0))l"^. (91) 


symmetry these expressions have to be modified. One 
can separate seven different isotopic projections 


Of{x;zi,Z 2 ) = {Pf)%(p°-'(x + zin)ip’’' (x + Z 2 n), (94) 


where / = 1,..., 7. The projectors Pf can be found in 

Ref. [itjl j. The one-loop anomalous dimension for the 
( f) 

operator O) is given by 



i ^ ^ 

e'v' {j + i){j + 2))^ 


(95) 


where A/ are the eigenvalues of the operator = 
daa'cfjbb'c invariant subspaces, RP/ = XjPf. The 

explicit expressions for A/ can be found in [l^. We note 
also that As = Ai and Ad = 2 A 3 . 

The modifications of the expressions (1711) . (I5T|) . (i5Bl) . 
(IC.13|1 and (1C.141) for the case of arbitrary projections, 
Pf, are the following: 


Here Zj is the one-loop renormalization constant for 
the operator Oj, 


u As 

7(j + 1)(j+2) 


(92) 


— Replace 7 ]^^ —>■ in all expressions. 

— Replace As —>■ A/ in the expressions for T^°\ 

rf yj6,A) rj,SE^ ^ 

— Replace As Ad ^ 2vf in the expression for 


and 

^1 = 1-^+0(^2), Z^ = l-^ + 0{u^). (93) 
12e 4e 

We have put the superscript (e) to the correlator in or¬ 
der to stress that even the tree-level correlator depends 
on e through a space-time dimension d = 6 — 2e. 


where izf are the eigenvalues of the invariant operator 
see Ref. 0. 


Representing the ratio of the correlators in the form 


t/m = ") + ...) 


(96) 


6 Results 


one obtains for the coefficient T; 


/( 2 ). 


The coefficient T) ’ in the ratio of the correlators is 
given by the sum of terms in Eqs. (|711) . (1511) . (I5B1) . (|C.13[) 
and (1C.141) . In the case of operators with other isotopic 




( 2 ) 


ab 


(97) 
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where 


r ( 2 ) = _1 (11 


= — ( - + 
^3,sd 12 " ^ 


2f + 7j-l 

72 V 3 ■ (j + l)(j + 2 )(j + 3) 
2f + 5j + l 


= 2 


3 ' (j + l)0- + 2)(j+3) 
f + 8f + lOj + 1 


jJf (^• + i)3(_^- + 2)3(j+3)’ 


A fri - 


1 


0 - + l)(j + 2 ) ^ 


7^(2) _ _ 

- 


and 


^( 2 ) _ _ 


3(j + l)(j + 2) 


(j + l)(j + 3) 

5^+3 - 4 

+ 2j + 5 

(j + l)(j + 2)(j + 3) 


3J (j + i)2(^- + 2)2 ■ 


(98) 


(99) 


Finally, comparing (I96p with the r.h.s. of Eq. (HU we 
get for the anomalous dimension 

7/( m *) = + ..., (100) 

where 7 )^ ^ ' is given by Eq. dMl) and ' has the form 


/( 2 ) ^ 27 ( 2 ) + -- 

^3 + (j- + 1)2(^- + 2)2 


2 A 




3(j + l)(j+2) 

A/Ad 

ij + 1 )(J + 2 ) 


Sj+2 - 4 


5j+2 - 3 ■ 


f + j - 1 

Hj + 1 )^ 0 ’ + 2 )^ 

2j + 3 


2 (j + l)(j + 2 )J 

J + IJ 


This expression completely agrees with the anomalous 
dimensions reconstructed from the two-loop evolution 
kernels [l^. We have also checked that the large j ex¬ 
pansion of the anomalous dimensions 7 ^( 11 *) in terms 
of the quadratic Casimir 


1 


J = ( i + 2 - e + -7/ {u*) ) ( j -t -1 - e + - 7 ) (m*) 


1 


2 


contains only even powers of 1/J 


27, 26 


7 Summary 


We have calculated two-loop anomalous dimensions of 
the leading-twist operators in the model using the 
approach proposed in Refs. Sim. Formally this me¬ 
thod allows one to gain one order in the perturbation 
theory. However, this advantage is illusory since one has 




Fig. 6 NLO diagrams for the correlator of the divergence of 
conformal operators (9C)j"^(2:)d(A^"^(0)). 


to calculate diagrams including finite parts instead of 
the pole terms in the standard approach. There is also 
no essential gain in a complexity of calculations. 

Nevertheless, both the contributing diagrams and 
the methods of calculation are quite different in the 
two approaches. Therefore the calculation of anomalous 
dimensions performed in this approach could provide 
an additional check of the results obtained within the 
standard approach. Of course, going to the next order is 
only possible with the advanced methods of computer 
algebra. 

Acknowledgements We are grateful to Vladimir Braun for 
useful discussions. This work has been supported by Deutsche 
Forschungsgemeinschaft through grant MO 1801/1-1. 


Appendix A: Discrete series representations of 
the stt(l, 1) group 

At the leading order the light-ray operator 0{x] Z\,Z 2 ) 
transforms according to the tensor product of discrete 
series representations of the group sm(1, 1). In this Ap¬ 
pendix we recall their basis properties. The discrete se¬ 
ries representation of the su(l, 1 ) group, in the stan¬ 
dard notations [l^ , is defined on the space of functions 
analytic inside the unit circle, \z\ < 1 , 

Dt{g-^)f{z) = {bz + a)-^^f{z'), (A.l) 


where z' = [az + b)l(bz a), and g = € su{l, 1 ) 

and s is an integer or half-integer. The generators of 
the group have the form 

S- = -dz, So = zdz + s, 5+ = z^dz + 2sz. (A.2) 

The invariant scalar product 0,0 

{Dt{9)f,Dt{g)'il3)s = {f,ib)s 
is uniquely defined (up to unitary equivalence) 

(/, V')® = J d^z g-siz) f{z) i}{z) , (A.3) 
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where 


A..(z) = ^(l-|zpr-^ 0 (l-|z| 


The space of analytic functions on the unit disk with 
the scalar product (IA.3I) is called a holomorphic Hilbert 
space, Z?+, see for a review Ref. [22|. The powers of z 
form an orthogonal basis, {ek{z) = z^, k = 0,1 ,...}, 
in this space. 


^k)s — ^kra I |^/c I Is — ^rak 


r{2s)k\ 


(A.4) 


'r{k + 2s) ■ 

The unit operator (the reproducing kernel) has the form 


ICs{z,w) ='^ek{z)ek{w)/\\e’"\\l = {1 - zw) (A.5) 
k=0 

and for an arbitrary function / € D'^ the following 
identity holds 


f{z) = J (fw fisiw)ICs{z,'w) f{w). (A.6) 

We note here that all formulae (EH) - ESI) have a 
perfect sense for any s > 1/2. 

Finally, we give a relation that turned out to be very 
useful in the calculations 


J (fwfJ.s+e{w)ICs{z,w)f{w) = 

= J cfiw fj.s{w)ICs-e{z,w) f{w) + 0{e'^). (A.7) 


It follows immediately from (IA.6I) if one replaces s —>■ 
s + e and expands it in e. 


Appendix B: One loop scalar product 

The one loop correction to the scalar product (15^ 
is determined by Eqs. (P|). To find the solution we note 
that the one loop kernel can be represented in the 
factorized formlj 

HW = 7=12^^12 = J2iJ' 21 , (B.l) 

where 

•^12 = {d2Z2l) ^ , J'12 = (^ 12 ^^ 1 ^ 12 ) ! 

^21 = {dizi 2 y^ , T 21 = (^2^524)”' ■ (B.2) 

® Let US remark that the evoluion kernel can be identi¬ 

fied with the sl(2)—invariant 7^—operator for a special value 
of spectral parameter, = TZs^ = i,B 2 = iy = The fac¬ 

torization of is a consequence of a factorization property 
of the 7?.—operator f^ . see also Ref. [23l . 


In order to obtain (1131 it is sufficient to notice that 
is nothing else as the inverse Casimir operator, 

(-^( 1 ))-! = -d,d2zl, = Ty7~l = 7yT~l. 

While the operator 1 —(g) 

the T operators intertwine the representations with dif¬ 
ferent spins. Namely, 

Ti2 ® Ti2 , 

yi2D^y (B.3) 

and similar for J^ 2 i- It results in the intertwining rela¬ 
tions for two-particle generators, = syysy\ 

V 21 (B.4) 

and so on. 

Next, we introduce one particle operator 
W^^y{z) = [f{z) - f{az)) , (B.5) 

such that = ('0{n + 2s) — '0(2s))z". This oper- 

(s) 

ator commutes with the generator 5'g while 

=-z. (B.6) 


Now let us check that 

+ Ku,iyi2w!y^yi2 + y2iw[^'>y2i) (b.7) 

gives solution to Eqs. (ESD- The first equation is ob¬ 
viously satisfied, [S'o°\ = 0. Next, making use of 

Eqs. dUl), dES]) one gets for the commutator 

= -{e - u^y^'^){zi +Z 2 ) 

— AgU* ^.Fi2Z2 J '12 + J21^iJ21^ ■ (B.8) 

Finally, one casts the r.h.s. into the necessary form tak¬ 
ing into account that [-J 12 , Z 2 ] = [-J 21 , zi] = 0. 

The scalar product can also be written in the form 




(^)- 


(J'l2/,lTi^4i2^)3 1 


+ (J'2i/,IE/^^J'2i4ia 


+ 0(e'), (B.9) 


where s* = 2 — e + y^s the conformal spin of the basic 
field at the critical point and (/, '!/')si.s 2 stays for the two 
particle scalar product. We have to mention here that 
the solution of Eqs. (|55)) is not unique. For instance, at 
one loop order w' = + Z, where Z is an invariant 

operator, [Z, 4°^] = 0, also satisfies Eqs. (I55|) . 
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Closing this section we give the standard represen¬ 
tation for the conformal operator, 




Z 1 —Z 2 —O 


(B.IO) 


The operator Oj is completely determined W a poly¬ 
nomial Z 2 ). It was known a long ago 23 that at 

the leading order 

P]{Z1,Z2) {Zi + Z2YCfl‘^'^ 

\Zl+ Z2 

C3 1‘2'\ 

where Cj is the Gegenbauer polynomial. Beyond the 
leading order one derives from Eqs. (IMl) and (Ell) 

Pj{zi,Z2) = {Zi + Z2y I P^^\zi,Z2) 


0 I 

-4u*As ^ hlp^^\zi,Z2)\ , (B.ll) 

fc=0,2... ^ 

where A = 2s* — 1/2, 
pf\zi,Z2) = 


j\r{ 2 \)r{\ +( zi — Z 2 


_ (jy _ 

r{j + 2X)r{j + x + ^) ^ \zi + z 2 j' 

and the expansion coefficients have the form 


(j + l)0- + 2) 


S' 


- 25 ,+1 - 


1 


{2k+ 3) 


(j + l)(j + 2) 

(fc + 1)! 


{j-'k){j + k + 3){jUy.{j + i)- 

The expression (IB.llIl agrees with the expression for 
the conformal operator obtained in 11 M- 


Appendix C: Loop diagrams 


We will present an answer for a diagram minus 
the counterterm AD^°''> in the form 

x'^inn) (£>(“) - = ul ■ KjTl°-\of\x)of'^)o . 

lf\x)0^^ 

diagrams in Fig. [S] and [5] we obtain 


Here {0j\x)0^^’)o is the tree level correlator. For the 


n{5,A) _ XsXd 


12 


^2j+2 


- Sj+3 - Sj+2 


T. 


(5.B) 


XgXd 


^ {j + l){j+2) 

(5,c) _ 2AsAd 


5 . 


2j+2 


-s. 


1 


'j+3 ' 


^ {j + l){j + 2) 

rp{b,D) ^ _ AsAd 


S2i+2 5j+3 -Sj 


3 + 1 . 
1 
2 


{j + 1)20- + 2f 


-1 , 
(C.13) 


and 


rji(6,A) _ 
3 




(1) J Ai) 


7 


2A,, 


52j-|-2 — 5j+3 — 5^+2 -I- - 


{j + l)(j + 2) L 


Sj+2 - 


T, 


(6.B) 


(i)-)2 J^ + 3j + 4 
2^^^' ^ (j + l)(j + 5)- 


{2j + l)(4j + 7) 

3(j + l)(j + 2) 

(C.14) 
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